Precision Measurement of the Df* — Mass Difference 

We have measured the vector-pseudoscalar mass splitting M(D* + ) — M(D^) = 144.22 ± 
0.47 ± 0.37 'MeV , significantly more precise than the previous world average. We minimize 
the systematic errors by also measuring the vector-pseudoscalar mass difference M(D*°) — 
M{D°) using the radiative decay D*° D°j, obtaining [M(D* S + ) - M(D+)] - [M(D*°) - 
M(D )] = 2.09 ± 0.47 ± 0.37MeV. This is then combined with our previous high-precision 
measurement of M(D*°) - M(D°), which used the decay D*° — > D°7r°. We also measure 
the mass difference M(D+) - M(D+) = 99.5 ± 0.6 ± 0.3 MeV, using the <f>ir + decay modes 
of the D s and D + mesons. 



I. INTRODUCTION 

Mass splittings between states with the same quark content but different spin configurations give essential 
information on the nature of the interquark potential. For example, masses of states with orbital angular 
momentum can be used to probe the contributions to the potential from spin-orbit and tensor forces between 
the quarks. High-precision measurements of mass splittings between states without orbital excitation (e.g., 
the pseudoscalar and vector S-states studied in this paper) give information on the relative contributions 
of chromoelectric and chromomagnetic terms to the Hamiltonian. Comparing the masses of resonant states 
having the same spin and charge configuration, but differing in the mass of one of the constituent quarks can 
isolate the effects of individual terms in the interquark potential. Of particular interest here is the vector- 
pseudoscalar mass splitting for cs mesons compared with cd mesons, as these are identical in the flavor SU(3) 
limit. Differences between them are presumably due to differences in the chromomagnetic contribution to 
the interquark potential and to the different value of the wave function at the origin, because of the different 
light quark mass. 

Using the decay modes D* + — ► D + ir° and D*° — ► D°it , CLEO II recently produced the definitive 
measurements of the mass splittings between the vector and pseudoscalar non-strange charmed mesons: 
M(D* + ) — M(D + ), and M(D*°) — M(D°) [?]. These high precision measurements were made possible by: 
a) the large data sample accumulated by the CLEO II experiment, b) the CLEO II crystal calorimeter, which 
allowed us to reconstruct the decay mode D* — ► Dir° with high efficiency and good resolution, and c) the 
fact that the decay through pion emission is close to threshold, giving excellent precision on the D* mass. 

Although there were a comparable number of observed events corresponding to the radiative decay D*° — ► 
D°j, this mode was not used because the larger Q-value degrades the mass-difference precision relative to 
the D*° — ► D°it mode. However, having measured the D*° — D° splitting to an accuracy of better than 100 
KeV using the pionic mode, we can use this to calibrate the mass difference measurement in the radiative 
mode. This, in turn, can be used to eliminate many systematic errors in our measurements of D* + — ► D s j. 

II. DETECTOR, DATA SAMPLE, AND EVENT SELECTION 

The CLEO II detector is a general purpose solenoidal magnet spectrometer and calorimeter. Elements of 
the detector, and performance characteristics, are described in detail elsewhere [?]. The detector is designed 
to have high efficiency for triggering and reconstruction of both leptonic and hadronic events. Charged 
particle momentum measurements are made with three nested coaxial drift chambers consisting of 6, 10, and 
51 layers, respectively. These chambers fill the volume r=3 cm to r=l m, where r is the radial coordinate 
relative to the beam (z) axis. Eleven of the layers in the main 51-layer drift chamber have sense-wires 
which are slanted relative to the beam axis to give measurements of the coordinate along z. More precise 
measurements of the z-coordinate are obtained from cathode pads located at the interfaces of the three 
tracking chambers. The system achieves a momentum resolution of (6p/p) 2 = (0.0015p) 2 + (0.005) 2 , where 
p is the momentum, measured in GeV/c. Pulse height measurements in the main 51-layer drift chamber 
provide dE/dx resolution of 6.5% for Bhabhas, giving good ir/K separation up to momenta of 700 MeV/c. 
Outside the central tracking chambers are plastic scintillation counters which are used as a fast element in 
the trigger system and also give particle identification information from time of flight measurements. The 
scintillation counters have a resolution of 154 ps as measured for hadrons, allowing better than 3<r ir/K 
separation up to momenta of 1.2 GeV/c. 
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Beyond the time of flight system is the electromagnetic calorimeter, consisting of 7800 thallium doped 
Csl crystals [?]. The crystal array gives an energy resolution of approximately 4% at 100 MeV and 1.2% at 
5 GeV. The central "barrel" region of the detector covers a solid angle of 75% of 4-7T. The endcap regions 
extend the solid angle coverage of the calorimeter to 95% of 4-7T, although with poorer energy resolution 
than the barrel region. The tracking system, time of flight counters, and calorimeter are installed in a 1.5 T 
superconducting coil. Flux return and tracking chambers used for muon detection are located immediately 
outside the coil and in the two endcap regions. 

A. Data Sample 

The data sample consists of 1.70/& -1 of e + e~ annihilations collected at CESR at energies just above and 
below the T(4S) resonance, and on the T(4S) resonance itself; the total data sample corresponds to about 
2 x 10 6 produced cc pairs. All events with 3 or more charged tracks, 1.5 GeV of energy measured in the 
calorimeter, and having a measured event vertex along the z-coordinate within 5 cm of the known interaction 
point, are accepted as hadronic event candidates. These events are then used for reconstruction of charmed 
mesons. 



B. Charged Particle and Neutral Particle Selection 

Our study requires the reconstruction of D s and D + in the <f>'ir + mode (with <f> — ► K + K~), and D° in 
the A'~7r + mode, as well as reconstruction of photons from the radiative transition between the vector and 
pseudoscalar states. We impose cuts on candidate tracks, requiring mainly that they come from the primary 
vertex. Candidate charged and neutral particles must satisfy the requirements listed in Table ??. We impose 
a 7T° veto on each photon candidate. The 7r° veto is implemented by matching photon candidates with other 
photon candidates passing the same quality cuts listed in Table ??. If their invariant mass falls within 2.5<r 
(approximately 12 MeV) of the known 7r° mass and if they give a good kinematic fit to the 7r° hypothesis, 
these photons are eliminated from further consideration. 

III. STUDY OF THE MASS DIFFERENCE 

We begin our study by focusing on the reconstruction of the c/>7r + decay mode. We require that the two 
candidate kaons from the <f> have particle identification information consistent with that expected for real 
kaons. There is a large background due to uncorrelated <f> and 7r + candidates, which peaks at cos6^ = 1, 
where 9^ is the decay angle of the <f> measured in the charmed-meson rest frame with respect to the charmed- 
meson momentum vector in the lab frame. The cut cos6^ < 0.8 is effective in reducing this background while 
retaining 90% of the isotropic signal. 

In decays of pseudoscalar charmed mesons into </>7r + , the <f> is polarized and its decay helicity angle (defined 
as the angle between one of the daughter kaons and the parent charmed-meson in the <f> frame) follows 
a cos 2 6 hgUcity distribution. To improve signal-to-noise we require that \cos6h e iicity\ >0.4. To suppress 
combinatoric background, we take advantage of the characteristic hard fragmentation function of charmed 
particles and impose the requirement x p (= "p candidate I Vmax) >0.5. With the above cuts, we obtain the c/>7r + 
invariant mass plot shown in Fig. 1; the mass plot has been fit to two Gaussian signals (representing 
D + — ► </>7r + , and D s — ► c/>7r + ) on top of a smooth background. The fit to the D + — ► 4>Tt + and D s — ► c/>7r + 
peaks yields approximately 400 and 1400 events, respectively. 

We use these fitted signals to determine the mass difference between the D s and D + mesons. Although 
uncertainties in the overall mass scale are on the order of 1-2 MeV, we expect the systematic error in 
the determination of the difference in masses to be much smaller. Contributions to the overall systematic 
uncertainty may arise from fitting (which we determine to be 0.25 MeV by varying the fit interval and 
the background function), and from possible differences between the lab momentum spectra of the <f> and 
7r daughters in the two cases. We probe the latter effect by fitting the mass difference in bins of scaled 
momentum x p , as shown in Fig. 2. The data are consistent with no variation as a function of x p at the 
67% confidence level, and we attribute a systematic error less than 0.1 MeV due to such a dependence. We 
arrive at a total systematic error of 0.3 MeV, and are therefore able to determine the difference in masses 
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FIG. 1. Invariant mass of (j>ir combinations. The smaller peak is the D —>■ 07r signal and the larger peak is the 
Df — > cf>w + signal. 
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FIG. 2. Mass difference between D and Dj , where both mesons are observed in the (j>ir mode, as a function of 
scaled momentum x p . 
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between the D s and D + to be 99.5±0.6±0.3 MeV. This value compares well with the present Particle Data 
Group value of 99.5±0.6 MeV [?]. 

IV. DETERMINATION OF THE D* s + -D s MASS DIFFERENCE 

With our large sample of Df — ► </>7r + , we can make a precise measurement of the D* + -D s mass difference. 
The previous best measurement of this mass difference was made by the ARGUS collaboration who obtained 
a value of 142.5±0.8±1.5 MeV [?]. The statistical precision of their measurement was limited by their small 
sample of D* + — ► D s j events in which the photons converted to e + e~ pairs. Their systematic precision was 
limited by the relatively large uncertainty in the calibration mode D*° — ► D°j. 

The CLEO measurement is made in the following manner. First, we combine Df candidates with photon 
candidates and use the resulting D* + — ► D s j signal to measure the mass-difference A]=M(D* + )-M(Df), 
where the s subscript on A indicates that we are considering the cs meson, and the 7 superscript indicates 
that the measurement is made using photon transitions. This raw mass-difference is still susceptible to errors 
in the overall photon energy calibration 1 which may be effectively eliminated as follows. Using the photon 
transition D*° D° 7, we similarly measure = M(D*°) - M(D°), which allows us to calculate the 
difference between the two mass-differences 8M=A^ — A^. This may then be used with the high precision 
measurement of AJ (using D*° -> _D°7r°) [?] to obtain A s = AJ + 6M. 

By imposing the same photon requirements in our measurements of the two radiative transitions under 
consideration we can extract A s relatively free of uncertainties in the absolute photon energy calibration. 
This technique is limited largely by differences in fitting the two signals due to the presence of the large 
D*° — ► D°it feed-down in the D°j mass-difference plot. There are no hadronic decays of the D* + states to 
produce such a reflection in the D* + — ► D s j mass-difference plot. 

A. Measurement of A] 

As discussed above, we reconstruct D* +, s in the mode D* + — ► D s j. In order to improve the signal-to-noise 
we cut on the decay angle # 7 of the photon in the D* + frame. Requiring cos# 7 > -0.7 eliminates a significant 
background, as is evident from Fig. 3. This requirement is made in addition to the other photon cuts 
detailed in Table ??. 

Figure 4 shows the distribution we obtain for M(4>Tt + ^) — M(4>Tt + ). The mass-difference distribution is 
fit to the sum of a smooth polynomial plus a "Crystal Ball Line Shape" 2 around the region of the expected 
signal. 3 The width of the signal and the magnitude of the tail are set at values obtained from Monte Carlo 
simulations. The mass difference we obtain from this direct measurement is 144.70±0.42 MeV, where the 
error is statistical only. 

B. Measurement of AZ and Determination of 6M 

We reconstruct D 0, s in the mode D° — ► K~ir + . The mass-difference signal M(D°j) — M(D°) is shown 
in Fig. 5 with two different fits overlaid. To obtain this mass-difference plot, we used the same photon cuts 
as in the D* + — ► D s j analysis. As before, we perform a fit (Fig. 5a) using the Crystal Ball Line shape 
function plus a smooth background. We explicitly exclude the low mass enhancement from D*° —> D°it 
from the fit region. The mass difference obtained is A^ = 142.61 ± 0.21 MeV (statistical errors only). This 



J The photon energy calibration is based on fitting the observed 7r° mass peak over a wide range of 7r° momenta. 

2 The Crystal Ball Line Shape is a nearly Gaussian distribution with a tail on the low end to take into account 
processes which may give an undermeasurement of the true photon energy. 

3 The enhancement at low mass-difference arises from misidentified D*-\ ► D + w° events where the D + decays to a 

three-body final state such as K~w + w + . When one of the final state particles is misidentified, kinematic reflections 
can occur in a mass region around the D s — > </>7r + signal. This has been verified by examining mass-differences using 
cf>w + combinations from the D s sideband region. 
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FIG. 3. M(07r + 7)-M(07r+) vs. cos6> 7 where Q 1 is the photon emission angle in the 07r + 7 frame relative to the 07r + 7 
direction in the lab. Transition photon candidates are required to have cos# 7 > -0.7, as described in the text. 
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FIG. 4. Mass difference between D s j and D s with fit overlaid. 
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may be compared with the value obtained from the tt° transition, A£ = 142.12 ± 0.05 ± 0.05 MeV, indicating 
that our overall photon energy calibration is understood to within 0.5% for the photons of interest in this 
measurement. To test fitting systematics, we perform an additional fit to this mass-difference plot (shown 
in Fig. 5b), where we explicitly account for the reflection from the hadronic mode. This is detailed further 
in our discussion of systematic errors. 

Comparing the mass differences obtained from Fig. 5a and Fig. 4, we determine 8M=A] — = 2.09±0.47 
MeV as summarized in Table ??. Combining this value with gives A s = 144.22±0.47 MeV. The errors 
quoted in both numbers are statistical errors only. 

C. Systematic Errors 

Systematic uncertainties arise from sources which affect A] and A^ differently and therefore introduce 
shifts in 6M. To the extent that the D* + and D*° fragmentation functions are different, photon energy 
calibration uncertainties can introduce systematic shifts, although the good agreement between A^ and AJ 
indicates that the photon energy scale is relatively well-understood. As is evident from Figs. 4 and 5, the 
background shapes are different in the two cases and there are therefore additional uncertainties arising from 
signal extraction systematics. 

We have studied possible biases using Monte Carlo simulations. Given input values of M(D*°) — M(D°) 
and M(D* + ) — M(Df ) we are able to recover values which are consistent with the input numbers after 
processing the Monte Carlo data through our analysis software. For the D* + — ► D s j transition, for example, 
inputting a mass difference between D* + and D s of 142.60 MeV, we recover a value of 142.55±0.15 MeV. 

We have investigated the dependence of the measured mass-difference on the photon energy and on the 
momentum of the D* + , which is correlated with the photon energy. Figure 6 demonstrates that the 
dependence of the measured mass-difference on transition photon energy is not large. Figure ?? shows the 
measured mass difference as a function of the scaled momentum x p of the D* + . The plot is consistent with 
no variation of mass difference with momentum. We therefore attribute no additional systematic error to 
such sources. 

There is also an uncertainty of ±0.5% in the absolute photon energy calibration which results in an error 
of ±0.7 MeV in A] and A^ as shown in Table II. However the contribution to 8M=A] — A^ is only ±0.02 
MeV since the systematic errors essentially cancel each other. 

Although systematics due to uncertainties in the overall energy calibration largely cancel, fitting system- 
atics remain. For the signal parameterization, we have checked that variations of signal shape produce shifts 
in both A] and A^ which track each other and therefore cancel in the value of 6M . The presence of the 
low mass enhancement due to the hadronic decay D*° — ► D°it can distort the shape of the background 4 
in the case of the calibration mode D*° — ► D°j. We have done a variety of fits using different assumptions 
for the photon line shape plus the possible background shapes in order to quantify the extent to which the 
hadronic decay can change the value of the mass difference we derive. Such a distribution is shown as the 
overlaid histogram in Fig. 5b. In this case, we have fit our observed signal to a sum of three pieces: a) 
a mass-difference background (whose shape is obtained from Monte Carlo studies) due to feed-down from 
D*° — ► D°it , 7T° — ► 77, where one of the 7r° daughter photons is reconstructed and the second is not detected 
in the calorimeter, b) a signal representing D*° — ► D°j, whose shape was also determined by Monte Carlo 
simulation, and c) a mass-difference background, obtained from M(K~ tt + j) — M(K~ tt + ), where the K~ir + 
combination is taken from the D° sideband regions. This gives a good fit to the data, indicating that we are 
able to account for the various components of the observed mass-difference plot. From this fit (Fig. 5b), we 
obtain a value of the mass difference A^ of 142.75±0.24 MeV. This compares well with the mass difference 
of 142.61±0.21 MeV obtained from Fig. 5a. We assign a systematic error contribution of 0.3 MeV to the 
measurement of A^ and 0.2 MeV to A], and conservatively assume the errors are totally uncorrelated in 
determining the contribution to the overall systematic error in 6M . 

The results of these measurements are summarized in Table ??. 



4 Note, however, that the hadronic mode is kinematically prohibited from producing background in the region of the 
D*° — > D°~f mass-difference signal. 
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FIG. 5. D*° - D° mass difference distribution, a) Shows a fit to the signal expected from true D*° — > D° 7 plus 
a smooth background, as done with D* + — > D s 7, and b) shows a fit to contributions arising from true D*° — > -D°7, 
D*° — > D°w°, and random photon plus fake D° combinations, as described in text. 
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FIG. 6. D*° — D° mass difference as a function of photon energy. 
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V. DETERMINATION OF UPPER LIMIT ON Z>*+ WIDTH 



It is straightforward to determine a limit on the intrinsic width of the D* + meson. The measured upper 
limit on the intrinsic width of the D*° is T < 2.1 MeV [?]. If we perform a free fit to the mass-difference 
signals observed in D*° — ► D°j (Fig. 5) and D* + — ► D s j (Fig. 4), using the same signal shape but allowing 
the width of the photon peak to vary, we obtain values of 4.50±0.24 and 4.29±0.40 MeV for the widths of 
the respective signals. Relating this to the intrinsic and experimental widths of the two resonances, we have: 

\J°d-<> + v 2 expU = 4.50 ± 0.24MeV, ^7 + ^ = 4.29 ± OAOMeV. (1) 

Assuming that the experimental resolutions u exv ti are identical for D*° — ► D°j and D* + — ► D s j, we can 
square and subtract these two expressions to obtain To* < 4.9 MeV at 90% confidence level. This technique 
is, at present, limited by the statistical precision on the A 7 measurements. 



VI. SUMMARY 



We have made a new measurement of the mass difference between the D s and the D + mesons, obtaining a 
value (99.5±0.6±0.3 MeV) in good agreement with the present world average, and with comparable errors. 
Calibrating our D* + -D s mass-difference using the mass-difference observed in the D*° — ► D°j mode, we 
determine 8M=A] — A 7 ,=2.09±0.47±0.37 MeV. Combining this value with our previous measurement of the 
D*° - D° mass difference [?], we determine M(_D s *+)-M(_D+)=144.22±0.47±0.37 MeV. This value is much 
more precise than the previous world average of 142.4±1.7 MeV [?]. 

It is of interest to compare the vector-pseudoscalar mass splitting for the cs system with that of the cd 
system. Two factors in the expression for the mass difference depend on the mass of the light quark: (i) the 
chromomagnetic effect is expected to be smaller for the cs system due to the heavier strange quark, but (ii) 
the square of the wave function overlap at the origin is expected to be larger because of the larger reduced 
mass of the strange quark. Our measurements indicate a larger vector-pseudoscalar splitting in the cs system 
than in the cd system, indicating that wave function overlap is the dominant effect. 

Finally, using the signal we observe in both the D*° — ► D°j and D* + — ► D s j modes, we determine the 
intrinsic full width of the D* + to be <4.91 MeV at 90% confidence level. 

Table ?? summarizes the vector and pseudoscalar splittings obtained by this and previous measurements. 
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FIG. 7. D*° — D° mass difference as a function of scaled D° momentum (x p ). 
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Measured D°(^ K~w + ) mass 
Charmed meson momentum 
|cos#| for photon candidates 
Photon candidates 
Photon shower isolation 
Photon energy 
Photon lateral energy deposition 
Photon 7T° veto 



<5 mm 
<5 cm 

< 2cr /3cr from expected for K/ir 
±2.5cr (10 MeV) of known mass 
±2.5(7 (21 MeV) of known mass 
±2.5(7 (26 MeV) of known mass 
x p > 0.5 
<0.7 (barrel region) 
unmatched to charged tracks 
>50 mrad from other showers 
> 50 MeV 

99% probability of coming from true photons 
inconsistent with coming from 7r° — > 77 



TABLE I. Summary of Cuts used in the Analysis, 
the interaction point.) 



("DOC A" denotes distance of closest approach to 



Raw AM 
Statistical Error 
Signal width systematic 
Signal tail systematic 
Momentum cut systematic 

cos8~/ cut systematic 
Background fit systematic 
Absolute _E 7 calibration systematic 



A7 
144.70 
±0.42 
0.1 
0.09 



0.3 
±0.7 



A 7 
142.61 
±0.21 
0.03 
0.06 



0.2 
±0.7 



8M-- 



=A7 - 

2.09 
±0.47 

0.1 
<0.05 



0.36 
±0.02 



A? 



TABLE II. Summary of Mass-Difference Results. (All numbers are in MeV.) 
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M(D+ - D u ) 
M(D+-D+) 
M(D*+ - D*°) 
M(D* S + -Dt)-M(D* - D°) 
M(D*+-D+) 



4.79±0.10 MeV [?,?] 
99.5±0.67 MeV [?],[this measurement] 
3.32±0.08±0.05 MeV [?] 
2.09±0.47±0.37 MeV [this measurement] 
144.22±0.47±0.37 MeV [this measurement] 



TABLE III. Summary of Charmed Meson Mass Splittings 
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